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Abstract 

We point out some similitudes between the statistics of high Reynolds num- 
ber turbulence and critical phenomena. An analogy is developed for two- 
dimensional decaying flows, in particular by studying the scaling properties 
of the two-point vorticity correlation function within a simple phenomeno- 
logical framework. The inverse of the Reynolds number is the analogue of 
the small parameter that separates the system from criticality. It is possible 
to introduce a set of three critical exponents; for the correlation length, the 
autocorrelation function and a so-called susceptibility, respectively. The expo- 
nents corresponding to the well-known enstrophy cascade theory of Kraichnan 
and Batchelor are, remarkably, the same as the Gaussian approximation ex- 
ponents for spin models. The limitations of the analogy, in particular the lack 

of universal scaling functions, are also discussed. 
PACS numbers: 47.27.Jv, 05.40.-a, 05.70.Jk 
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I. INTRODUCTION 



Turbulence is characterized in a generic way by self-similar spectra, a feature that is 
also found in other physical situations, such as critical phenomena or random walks. In real 
space, this property is usually the signature of well-known statistical mechanisms such as the 
inverse power law decay of two-point correlation functions and the non-Gaussian behavior 
of the relative dispersion of random variables of the problem. Some attempts have been 
made to relate scaling in turbulence with critical phenomena and field theories More 
recently, experiments have shown that the power consumption in closed turbulent flows 
and the total magnetization of spin systems at criticality have similar, and to some extent, 
universal, probability distribution functions 0. The two-point velocity correlation function 
has been a quantity extensively studied in hydrodynamic turbulence, mostly under the 
related form of velocity increments 0. In practice, real space experimental and numerical 
analysis of eventual scaling laws is a difficult, space and time consuming task, because the 
asymptotic regimes are reached very slowly as the Reynolds number goes to infinity. On 
the other hand, relatively few studies have been devoted to the characterization of the 
spatial structure of the vorticity (cJ = V x u). Although, in Fourier space, the enstrophy 
spectrum of an incompressible fiow is very simply related to the kinetic energy spectrum, 
the vortex statistics and correlations in real space are not so simply related to their velocity 
counterparts, but have distinct properties. This difference is illustrated in a spectacular way 
in two-dimensional turbulence, where the enstrophy (cJ^(x)) cascades toward small scales 
while the kinetic energy (u^(x)) follows an inverse cascade to large scales. This is essentially 
due to the absence of vortex stretching in two dimensions. Because of this property, the 
enstrophy obeys a conservation equation, and (cJ^(x)) is dissipated even in the very low 
viscosity limit, contrary to the energy 0,^. 

A vortex is one of the simplest kind of structure generated by hydrodynamic instabilities; 
an isolated structure is likely to produce slowly decaying flow flelds in space (typically with 
1/r, dependence) and consequently long range velocity correlations. Physically, a flow 
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can be structured even if the vorticity vanishes nearly everywhere, as in the case of an 
assembly of point vortices in two dimensions. Hence, the study of vorticity distribution 
functions should reveal a deeper, less intuitive degree of order in flows. This occurs in 
two-dimensional turbulence, where the vorticity tends to concentrate in spatially extended 
structures, of particle-like character. Indeed, numerical studies of two-dimensional decaying 
turbulence have shown that the long time dynamics is dominated by a few large coherent 
axisymmetric vortices that contain most of the enstrophy DH. Dynamically, it is commonly 
believed that these structures are the result of many vortex mergers. Such flows can be 
often well described by a modified Hamiltonian, deterministic model of point vortices [§. 
These rather well isolated structures of well defined shape and size have the advantage of 
simplifying the picture of the fiow, because the organization of initially numerous incoherent 
vortices into blobs reduces the number of relevant degrees of freedom. 

In this study, following a route similar to the analysis by Babiano et al. of velocity 
structure functions, we characterize the vorticity structure in real space for two-dimensional 
homogeneous unforced turbulence, keeping in mind its possible representation in terms of 
assemblies of vortices. We consider the normalized vorticity autocorrelation function, and 
calculate it assuming a phenomenological self-similar energy spectrum in the inertial range, 
E{q) ~ q^^. In experiments or numerical studies, /i is usually found to be around 3 P], the 
value predicted by the Kraichnan |^ and Batchelor reference theory, hereafter referred 
to as KB, or steeper (/i ~ 4). 

In the next section, we show that the vorticity autocorrelation function exhibits distinct 
behaviors depending on the value of the exponent fi. In section III, we show that it is 
possible to draw an analogy between two-dimensional turbulence and systems of spins near 
criticality for values of fx larger than 1. The correlations can indeed be written under the 
form of a scaling law involving a diverging correlation length when the Reynolds number 
goes to infinity. The critical exponents that correspond to the KB theory precisely coincide 
with the Gaussian exponents of spin models or random walks. Conclusions including a brief 
discussion of the three-dimensional case are presented in section IV. 



II. VORTICITY POWER SPECTRA 



Let us consider the vorticity autocorrelation function 

(J(x) ■tJ(x + r)) 

that has been normahzed to one at r = 0. R^^ can be expressed as a function of the Fourier 
transform of the vorticity, 



tu(q) = J ^{^) e-'"^-"". (2.2) 

through the relation 

/rfq (cJ(q)-cJ(-q)) e^^- 

r = — — ^ — . 2.3 

/rfq (w(q) ■ w(-q)) 

The vorticity power spectrum (cJ(q) ■ a;(— q)) can easily be related to the kinetic energy 
spectrum E{q) of the flow u. In two dimensions one gets, in the homogeneous isotropic 



case, 



E{q) = ^ (u(q) • u(-q)) , (2.4) 

where V is the volume of the flow. If the flow is incompressible, one has the identity 
(d;(q) ■ u;(-q)) = g^(u(q) ■ u(-q)), so that 

{u{q)-u{-q))=47rV qE{q) . (2.5) 



Replacing the identity (|2.5| ) in expression ( p.3|) and integrating over the angular variable 
gives: 

R .. Idqq'Eiq)Mqr) 

^ !dqq^E{q) ' ^^'^^ 

where Jo is the Bessel function of the first kind. In the following, one assumes that the 
energy spectrum is self similar in the inertial range, which extends between an integral scale 
/ at which energy is (initially) injected and an enstrophy dissipation scale a: 
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E{q) ~ q-'', 27i/l <q< 2n/a . (2.7) 

The spectrum is supposed to vanish at length scales smaller than a. If the observation scale 
r is much smaller than /, one can assume in addition that E{q) = for q < 27r//. (We 
discuss in section IV the modifications introduced by considering the part of the spectrum 
that extends beyond the inertial range.) Note that, in the general case, the scale I, and hence 
the integral Reynolds number, may increase with time, as in freely decaying processes [Q. 
This is not a restriction to our analysis, which deals with instantaneous spectra. Relation 
( p.6[ ) simply turns into: 

RUr) = (£" dq q-^'^-'^Mqr)^ (£" dq g-^^-^)) ' , (2.8) 

where 

qi = 27r// 

(2.9) 

qa = 27r/ a . 

Note that the normalized velocity autocorrelation function, Ru{r), has the same expression as 
Eq.( p.8|) , replacing /i by {fx + 2). It can be seen that both the numerator and the denominator 
of Eq.( p.8|) diverge when / — oo if > 3. Hence, we expect the behavior of Ruj{r) for > 3 
to differ from its behavior for fi < 3, while the analytical form of Ru does not undergo an 
abrupt change at this point. In the following we discuss the different cases encountered. 

A. Case fi > 3 

To study the correlations in the interval r ^ a, we can replace the upper limits of the 
integrals of Eq. (|2.8|) by infinity. Thus, R^jij-) scales as a function of r// only, and can be 
rewritten as 

/•oo 

R^r) = 1 - ifi - 3)qr' dq q'^^'^^ [1 - Mqr)] , a <^ r. (2.10) 

"'<?( 

If 3 < ;U < 5, the integral of Eq.( |2.1(J| ) is finite as qi goes to zero. After the variable change 
X = qr, one obtains the first order expansion in r/l, 
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] , a<r</;3</i<5, (2.11) 



where F denotes the Gamma function. If 5 < /i < 7, we use the second order expansion 
of the Bessel function Jo{x) = 1 — /4 + O(x^), and replace in Eq.(|2lO|) [1 - Jo(x)] by 
[1 - x^/A - Jo{x)] + xV4. We obtain 

RUr)c:.l-^-^^{qir)'-C,{qirr-^ 5 < < 7 , (2.12) 

with 

POO 

C^={fx-3) / dx x-^^-^^ll - x^A - Jo{x)] . (2.13) 



The leading behavior of Rui{t) is given by the r^-term. It is easy to show that this property 
stays valid for all /i > 5: 

i?^(r)-l- f ~ (qirf , a < r < /. (2.14) 
4(/i- 5) 



The results displayed in Eqs. ( p.ll|) and (|2.14| ) are the analogues of those obtained by Babiano 



et al. 1^ for the velocity structure functions (in the cases 1 < /i < 3 and /i > 3 respectively). 
The power-law exponent of the corrections for small separation distances does not depend 
on /i for /i > 5. The first order expansions ( p^.ll[ ) and ( ^.141 ) show that the characteristic 



length ruling the decay of R^j is ^. We conclude that the vorticity correlation length is of 
order of the integral scale /. 



B. Case 3/2 < < 3 

For energy spectra less steep than —3, the denominator of expression (|2.8|) is well defined 
as / goes to infinity, but diverges when a goes to zero. In turn, if 3/2 < /i < 3, the numerator 
of Eg. (p78|) does not necessitate infra-red nor ultra-violet cut-offs to be finite, provided that 
r > 0. This implies that in the limit where I /a is large, Ruj{t) can be much smaller than 
one even if r <C /. 
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If a ^ r ^ /, one can set qi = and ga = oo in the numerator of Eg. ( p.8|) . Using 



properties of Bessel functions [11|, one obtains: 



R.{r) ^ - ) . (2.15) 



a 



The correlation function decays algebraically with the distance r. The value /i = 3 separates 
two different scaling regimes of R^^. At the transition, the numerator of Eg. ( |2.8|) behaves 
like ln(//a), and making the variable change x = qa, one obtains after a simple calculation: 



T 

1 - Joix 



fi = 3. (2.16) 



When the separation distance is much larger than the dissipation scale, the correlations 
decay logarithmically: 

, , \n(l/r) - luTT - 0.577 

R^r) ~ \ ni ^ ' a < r < / ; /i = 3. 2.17 

ln(//a) 



C. Case n < 3/2 

For values of fi lower than 3/2, both the denominator and the numerator of Eg.( p.8|) 
diverge when a goes to zero. In order to avoid convergence problems due to oscillations 
when integrating the Bessel function Jq over a large but finite interval, it is now necessary 
to replace the discontinuous energy spectrum (|2.7|) by a continuous spectrum that models 



better the effects of viscous dissipation [12|, for instance: 



qr M g 9 /Qa a q > q, 
E{q) ~ <( ' (2.18) 
if g < qi. 

In Eg. (|2.18| ), the dissipative spectrum at large wavenumbers g > is described by a Gaus- 
sian function |T3[. Setting qi = 0, relation ( p. 61 ) yields: 

O POO 

Ru{r) = ^--^g;(3-A') / dq e-^'/^'q'-'^Mqr) , (2.19) 

r(^l Jo 

that can be recast as 



RUr)= , (2.20) 

where iFi is a hypergeometric function. For g^r ^ 1, and if /i 7^ 1, behaves as [|ll|] 

Mr)^-7—v[-) ■ (2.21) 



r(¥) 



V a 



If yU > 1, the above expression is positive and does not differ from the behavior displayed by 
Eg. ( p. 151) . In turn, if /i = 1, the prefactor of the algebraic dependence in Eq.( |2.21| ) vanishes 



and relation ( |2.2CI| ) simply reduces to a short range function: 

i?^(r) = e-(""/'^)'. (2.22) 
The energy spectrum /i = 1 is known to correspond to a flow composed by statistically 



independent point vortices [|T4|. In that limit, it is easy to show that the vorticity correlations 
reduce to 

1 if r = 

RM)\,=i = \ (2.23) 
if r > 



Replacing a by zero in Eg. ( p.22D , we recover the distribution given by Eq. (|2.23|) . 



If fi < 1, vortex correlations for large separation distances r are still given by relation 
( p.21| ). However, the prefactor T[{^ — l)/2]^^ is now negative. This case is physically very 



different from the situation /i > 1, where Ruj{t) is always a positive, slowly decaying function. 
Indeed, for < 1, Ruj{r) decreases for small r and has its first zero for r ~ a: this is the 
signature of a short range function, such as Eq.( p.22|) . We conclude that R^ is short range 



and characterized by a correlation length of order a in the range /U < 1. 

III. DISCUSSION 

The results derived in the previous section show that the two-point vorticity autocorre- 
lation function can be recast under scaling forms, such as 

R.ir) = ([y' f{r/0 ■ (3.1) 



The above relation defines an exponent t] and a correlation length ^. Pair correlations in 
systems near a critical point satisfy similar scaling laws, where the scaling function / decays, 
say, exponentially, and C, diverges as the temperature approaches the critical temperature. 
In the present context, the vorticity has various spatial structures depending on the slope of 
the energy spectrum. The correlation length ^ is small and equals the dissipation length a 
when /i < 1, see Eq.( p.22[) . The case /i > 1 is more interesting because, then, the correlation 
length is large: ^ = I ^ a, see Eqs. (|2.11|) and (|2.14|) . (In the range 1 < /i < 3, if the 
separation distance r is no longer small compared with /, it is easy to check that expression 
( p.l5| ) or ( p.21|) reduces to the form ( p.l|) with ^ = I.) These two distinct regimes are indeed 
driven by the behavior of the numerator of Eq. (|2.8|) with respect to the ultra-violet cut-off. 
When this cut-off is necessary in order to get a finite integral, the correlations are short 
range; otherwise, they are long range. This abrupt change is reminiscent of a similar feature 
encountered in problems of elastic interfaces at thermal equilibrium, namely the transition 



from flat to rough surfaces as the surface dimension decreases across the value d = 2 [15 



In that problem, the r.m.s of the height difference between two well-separated points can 
be seen as a correlation length. In high dimensions (analogous here to low values of fi), 
this difference is constant, of order of a small cut-off scale identifled with the thickness 
of the interface (here, the dissipation scale). However, if d < 2, the height difference is 
"macroscopic": it depends on the size of the system, i.e. the separation distance between 
the two points. 

The main scaling properties are summarized on the diagram of Figure 0. Flows with short 
range order can be pictured for instance as a gas, or liquid, of nearly point-like vortices, with 
radii of order of the smallest characteristic scale, the dissipation scale a. Using the language 
of spin models, these structures are also analogous to high temperature states, where a 
would correspond to the lattice spacing. The part of the diagram where /i is larger than 
1 corresponds to vorticity distributions with long range spatial order. In the following, we 
focus on this region and develop an analogy with critical phenomena. 
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A. Critical exponents 



An analogy between a physical problem and critical phenomena requires that some of its 
statistical properties take a simple scaling form described with a few relevant parameters in 
a particular asymptotic limit [0. This can be done, for instance, in the context of random 
walks, where the analogue of the spin-spin correlation is the probability of presence (or the 
end-to-end distance distribution), and the analogue of the temperature difference from the 
critical temperature is the inverse of the time (or length of the walk) |T7| , P^ . 

In turbulence, the integral Reynolds number of a flow is defined as Re = vI/uq, where 
V is the r.m.s. of the velocity and uq the kinematic viscosity. The ratio {l/aY, where d 
is the space dimension, is often seen as the number of degrees of freedom of the flow 0. 
Phenomenological theories involving dimensional arguments P,p!0| usually lead to power law 
behaviors oil/ a with Re, say 

Ir^aRe" . (3.2) 

Since / is the correlation length in the region > 1 (see Figure P, v can be understood 
as a correlation critical exponent. Indeed, as the Reynolds number goes to infinity, the 
correlation length ^ = / diverges in units of a, and the small parameter that separates the 
system from criticality is the inverse of the Reynolds number. 

The correlation function exponent r] of the scaling form (|3.1| ) defines a second critical 
exponent, = 3 — /i in the range 1 < /i < 3, and = for /i > 3. The transition (/i = 3) 
where correlations no longer decay algebraically at infinite Reynolds number coincides with 
the spectrum of the KB theory. At this particular point, Ri^ decays logarithmically as shown 
by Eg. ( p. 171) . Table (|) shows the main connections with critical phenomena, for instance 



with a system of Ising's spins at temperature slightly above T^. However, the correspondences 
are purely statistical: considering few body problems, it is obvious that vortex interaction 
and dynamics [|1^ can not be compared with spin interaction. Concerning this last point. 



the comparison with spin systems is much more precise in the case of a superfluid |^ 
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Different theories of turbulence will in general lead to different values of the exponents 

V and ri. In phenomenological theories, conservation laws and other arguments (like di- 
mensional arguments) determine the exponent of the energy spectrum, as well as prefactors 
depending on dissipation rates. Since the size of the inertial range must be consistent with 
the dissipation rates, once the shape of the spectrum is known and the energy injection scale 
is fixed, the dissipation scale a (and consequently the value of the exponent v in relation 
( p.2|) ) is unique. Hence, the critical exponents ri and v are generally not independent. 

The KB approach to two-dimensional turbulence assumes that the enstrophy flux through 
wave number q is independent of q at large q\ a dimensional argument leads to /i = 3 and 
{l/aY ~ [SI- The critical exponents associated to the KB theory are hence 77 = and 

V = 1/2. These exponents are, remarkably, the same as the exponents of the Gaussian 
approximation in spin models, or the Brownian motion exponents (see ll^). Although 
the meaning of this coincidence remains unclear, first, it can be noted that a mean field 
approximation in the present fluid dynamics context corresponds to the fact of neglecting 
the fluctuations of the enstrophy dissipation rate during the enstrophy cascade. Secondly, 
the KB theory is based on two conservation equations: for the kinetic energy (u^) and for the 
enstrophy {00^) ■ Such a description is incomplete: indeed, in two dimensions the vorticity 
follows the fluid motion in the inviscid limit, hence any continuous functional {F{uj)) is 
invariant with time, or must satisfy a conservation equation in the slightly dissipative case. 
Assuming the Gaussian approximation, (cu^"^^) = and (cu^") = (cu^)" for any integer n, the 
conservation of (cu^) implies the conservation of any functional F{uj) = J2n dn^^- Hence, the 
KB theory ressembles a Gaussian approximation, by focusing only on the second moment 
of the vorticity. 

B. Lack of universality 

Although microscopic details (here the dissipation processes at g > g^) disappear in 
the behavior of Ruj{t) for large separation distances, the analogy with critical phenomena 
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presented in this section encounters a limitation when looking more closely at the scaling 
function / in Eq.( |3.1| ). In the case of turbulence, / does not rigorously converge in distinct 
"universality classes", since it does not depend on the self similar energy spectrum of the 
inertial range only. The scaling function / essentially depends on the shape of the energy 
spectrum at low wave numbers. In general, the integral scale is not the largest scale in 
turbulent flows and coherent structures may emerge from external constraints: the wave 
numbers lower than 2ti/1 contain a non- negligible part of the energy and this part of the 
spectrum is driven by "non-universal" conditions like interactions between large eddies, 
depending on the system size and boundary conditions. Hence, only short distance scaling 
is generic in turbulence, as already noted by Eyink et al. [|l|. 
Let us consider the modified spectrum: 

^fe) J (3.3) 
yq ^ qi <q <qa- 

In two dimensions, the exponent fiQ = 1 corresponds to an equipartition kinetic energy 

spectrum at low wave numbers; the value /xq = 3 has also been proposed for decaying 

turbulence [0]. It can be easily shown that in the range r -C /, Eqs. (|2.15|) and (p.20|) 



remain unchanged with the spectrum given by Eq. (|3.3| ). However, if one wishes to compute 
the scaling function, in the case /i > 3 for instance, Eq. (|2.1(]| ) turns into 



RUr) = / (y) = 1 - [-j^ + ^) 9"^'"'^tl - "^o^^^)] ' (3.4) 

r < /; > 3. 

Hence, the differences introduced by the modified spectrum in the scaling function /(m), 
defined by Eq.( |3.1| ), already appear at the first order expansion in u: the prefactor is also a 
function of yUo- As u grows, the decay of f{u) depends more strongly on the "non-universal" 
part. 

However, with E{q) given by Eq.( |3.3| ), the leading behavior of f{u) for u ^ 1 is gener- 
ically exponential. A famous example is provided by the Ornstein-Zernike (OZ) approxi- 
mation for critical phenomena, where the spin-spin correlation function reads {S{Q)S{r)) oc 
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/ d^k exp(ik.r)/(/c^ + ^ ^) ~ r ^/^exp(— r/^) for r ^ ,^ corresponding in the present 



description to /ig = — 1 and /i = 3. Recent experimental results on two-dimensional decay- 
ing turbulence in soap films have shown clear evidence of spectra with both power-law and 
inverse power-law regimes, with respective slopes close to the values /xq = 1-5 and /i = 3 [If. 
Analytical calculations can be performed in the equipartition case /io = 1 and ^ = 3: let us 
rewrite the energy spectrum under the form E{q) ~ q/{q^ + C,~^) with ^ = I. The vorticity 
correlations read: 

R^{r) oc / dqq . (3.5) 

Jo + 4 

Differing from the OZ result, the asymptotic limit of Eq.( p.5D exhibits oscillations ||23 



f{u) oc exp(-M/v^) cos(m/ V2 + 7r/8) , m > 1. (3.6) 



We now outline a physical interpretation of the scaling function /. As mentioned be- 
fore, some two-dimensional decaying turbulent flows with a spectrum given by Eq. (|3.3|) are 
well pictured by an assembly of well separated vortex blobs with particle-like character. 
This has been, for instance, very clearly observed in some recent numerical simulations on 
homogeneous Rossby wave turbulence |Q. Let us assume that these vortices are all of 
size ^ and with the same circulation modulus. The vorticity distribution inside one blob 
is thus described by the short distance behavior of the vorticity correlation function, see 
e.g. Eq.( p.l7| ). In turn, the long distance behavior of R^^ accounts for the weak correlations 
between different vortices. Following a route similar to the one presented in ref. pS), it 



is easy to show from the definition (|2.1| ) that, if there are the same number of blobs with 
positive and negative circulation so that the flow has no net overall rotation, 

R.{r) ^^[gi{r) - g^{r)] , r > ^. (3.7) 

Using the notations employed for ionic liquids, gi represents the pair correlation function 
between two vortices of same circulation ("like"), and gu the pair correlation function be- 
tween two vortices with opposite circulation ("unlike"). The oscillating decaying shape of 
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particle pair correlation functions is a generic feature of liquids, in particular, those that are 
ionic [^]. The oscillations of i?^ obtained in a particular case through Eg. (p^) , although 
having a very different physical origin, qualitatively support the idea that the structure of 
homogeneous assemblies of coherent vortices can be "liquid" -like. 



C. Susceptibility 

Since the autocorrelation function is an averaged quantity, it cannot of course provide 
by itself all the details on the complexity of vortex structures observed in turbulence, and 
contained in higher moments. Yet, large coherent vortices can be seen as a particular 
manifestation of long range order up to the integral scale, similarly to magnetization domains 
in systems of spins near a critical point. These vortices produce large fluctuations of vorticity 
on scales smaller than the integral scale, and possibly spatial intermittency. The fluctuations 
of the total vorticity contained inside an observation window of area A are non Gaussian 
(for /i > 1) when A is lower than the square of the integral scale. To see this, let us define 
the circulation, 

^A = ^rJyj^Jj^-^{^) . (3.8) 

The standard deviation ^ {^"a) of Qa from its mean value {{^a) = 0) can be estimated by 
considering Qa as a sum of random correlated variables. According to the definitions ( p.l|) 
and (|3.8| ), one gets 

1/2 



(a JjrRUr)^ . (3.9) 
When is short range, one obtains the Gaussian behavior. 



(nl) ~ /i < 1, (3.10) 

at any scale larger than the dissipation scale. In turn, when the correlation length equals the 
integral scale /, the scaling exponent is non-Gaussian (i.e larger than 1/2) at scales smaller 
than I. With the help of Eqs. (lOI) and (|3.9|) , we indeed find: 
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{A<f) (3.11) 
with 

(u + l)/4 1 < /i < 3 

a = l-r]/A= \ (3.12) 

I 1 /i > 3. 

To derive a critical relation analogous to the one for the susceptibility, we invoke the argu- 



ments of Bouchaud and Georges [|T8[ in their statistical interpretation of the general relations 
between critical exponents. If the area A of the observation window exceeds one can re- 



place the integration domain of Eq. (|3.9| ) by multiplying the factor A by I'^/P, one gets: 

= (A//2)l/2^(]2^)l/2 _ (^/;2)l/2(;2)a ^ ^^^2 (g^^g) 

Hence, the system can be pictured as A/l"^ independent regions of size / at criticality. We 
define the susceptibility per unit 

x = ^. (3.14) 

When A is larger than the integral scale squared, x becomes independent of A. Taking into 
account the dissipation scale a in the calculations, we find from expression ( p.l3|) combined 
with Eqs.(|3) and (|3l2D : 



X~a2 , (3.15) 

with the exponent 7 satisfying the same relation as in critical phenomena: 

7 = z/(2-r7). (3.16) 

Relation ( p.l5| ) shows that the susceptibility diverges when compared with the "microscopic" 
susceptibility a^. Note that contrary to the lattice spacing in spin models, the dissipation 
length a may not remain constant while other parameters of the flow are varied. In any 
case. Re, (cu^) and a are independent variables. 
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IV. CONCLUSION 



We have studied the two-point vorticity correlation function in two-dimensional turbulent 
flows characterized by a self-similar energy spectrum within an inertial range limited by two 
very different length scales. If the slope is steeper than 1, the vorticity distribution presents 
some similitudes with the magnetization in models of spins near critical point. 

On the grounds of phenomenological theories of turbulence, the analogy with critical 
phenomena is made possible because the ratio between the integral scale (the macroscopic 
scale at which energy is injected) and the dissipation scale diverges as a power law of the 
Reynolds number in the large Reynolds number limit. The divergence of susceptibility 
is associated here with the increasing fluctuations of the circulation taken around a large 
contour, as the Reynolds number increases. As a surprising result, the (non-independent) 
critical exponents corresponding to the Kraichnan ^ and Batchelor JlOl (KB) theory are the 
same as those of the Gaussian approximation result in spin models. A possible connection 
between these two phenomena could be as follows: The KB theory only considers the first 
non- vanishing moment of the vorticity {i.e. {uj"^)) as an inviscid invariant of the motion, 
among all the infinite number of possible invariants associated with the vorticity. If one 
invokes the Gaussian approximation (u;^") = (u;^)" for any integer n, the conservation of 
(cu^) is indeed sufficient for any C°° functional F{uj) to be conserved. Another particularity 
of the KB theory is that the vorticity correlations decay logarithmically for distances larger 
than the dissipation scale and shorter than the integral scale. The slope of the KB energy 
spectrum (—3) indeed corresponds to a transition value, where the correlation function no 
longer decays algebraically at infinite Reynolds number. 

Flows with the same inertial range spectrum may have distinct scaling functions, de- 
pending on their large-scale configuration; the non-universality of correlations can already 
be noticed by looking at the prefactor of a first order expansion in r/^. However, these func- 
tions commonly decay exponentially at distances much larger than the correlation length. 
On the basis of experimental results , we argue that pair correlations between large-scale 
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coherent vortices may be liquid-like. 

Most of the statistical remarks made above are not qualitatively modified in three di- 
mensions. It can be shown that for a Kolmogorov —5/3 law, the vorticity correlations decay 
algebraically and satisfy relation ( |3.1D with 77 = 4/3 |^. In the general case, at finite 
Reynolds numbers, the circulation modulus introduced in relation ( ^.81 ) is a sum of random 
variables if the integrated area extends over many integral scales, see Eq.( p.l3|) . However, it 



scales anomalously in the inertial range, see Eq. (|3.11|) . This is probably a feature satisfied 



by other global quantities in turbulent fiows, like the power consumption studied experimen- 



tally in refs. P,|28[|. The present analysis is consistent with some of the results presented in 
these two references: in a closed fiow confined to a size of order of one energy injection scale, 
the power consumption fiuctuates strongly. In turn, in an open fiow, many large eddies can 
develop and the fiuctuations are Gaussian, according to the central limit theorem. 
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FIGURES 

FIG. 1. Scaling forms of the vorticity autocorrelation function in the space of the parameter 

/X. If the slope of the energy spectrum is lower than 1, R^^ is a short range function. In turn, wlien 
H > 1, the correlation length is of order of the integral scale I. Note that for /x > 3, the critical 
exponent r] is identically zero. 
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TABLES 



TABLE I. Two-dimensional turbulence and critical phenomena compared. 



2d turbulence 



Critical phenomena 



i?,(r) = (d;(x).d;(x + r))/(a;2(x)) 

1/Re 0) 

RUr) - {r/a)-y{r/l) 

Dissipation scale a 

Integral length: I ~ a{l/Re)~'^ 

X = A-^ ((/^ dx u;(x))^) ~ a\\lRe)-^ 

Exponents relation: 7 = i/(2 — 77) 

Kraichnan h Batchelor's theories: 77 = 0, z^ = l/2 



Spin correlations (S'(x)S'(x + r)) 
t= IT-Tel/Tc H 0) 
(5(x)5(x + r)) ~ r2-<i-V(r/0 
Lattice spacing a 
Correlation length ^ ~ 
Susceptibility: x ^ ((Z^i Sj)^^ ~ 
7 = 1/(2 - 77) 

Gaussian result: 77 = 0, z/ = l/2 
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R„(r) ~ r-" f(r/^) 



gas or liquid of 
small vortices 

^ ~ a 



long range order 



^ - 1 
?7 = 3-/z 



^ - 1 
r] = 



Fig. 1 



